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J> . Abstract 
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An analytical theory of low frequency electromagnetic waves in metallic photonic 
, crystals with a small volume fraction of a metal is presented. The evidence of the 

existence of such waves has been found recently via experiments and computations. 
We have obtained an exact dispersion equation for oj(k) and studied the cutoff fre- 
quency ujq = uj(0) as a function of parameters of the photonic crystal. An analytical 
expression for the permittivity e is calculated. It is shown, that if the crystal is 
embedded into a medium with negative /J,, it has no propagating modes at any fre- 
. quency. Thus, such a compound system is not a left-handed material (LHM). The 

£Z] ' recent experimental results on the LHM are discussed. 

o 
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Veselago has shown [1] that if in some frequency range both permittivity e and 
permeability \i are negative, the electromagnetic waves (EMW's) propagate 
but they have some peculiar properties, from which the primary one is that 
vectors k, E, H form rather a left-handed than a right-handed triple of vectors. 
Materials, with this property are named left-handed materials (LHM). 

The San Diego group [2] uses metallic photonic crystals (MPC's) as a techno- 
logical base for the LHM in the GHz frequency range. This idea comes from the 
pioneering computational and experimental studies of a few groups [3-7] which 
have found very low frequency EMW's in the MPC. These EMW's propagate 



1 We are grateful to A. Zakhidov and V. Vardeny who have attracted our attention 
to this interesting problem and to L. P. Pitaevskii and B. L. Spivak for important 
discussions. This work supported by the NSF grant DMR-0102964. 
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above a very low cutoff frequency for which various groups obtained differ- 
ent values. Since the parameters of the MPC's were also different, it was not 
clear whether or not the physics of these modes is the same. The suggested 
interpretations are very controversial. The group of Soukoulis[5] qualitatively 
interpreted the effect of propagation in terms of waveguide modes, while the 
group of Pendry[6] presented a completely original physical picture based upon 
a new longitudinal mode, called "plasma mode" . According to Pendry et al. [6] 
the resulting effective permittivity has a plasma-like behavior 



Lu(LU + iT) , 



however, expression for the "plasma" frequency contains the velocity of light 
and has a form 

2 _ 2ttc 2 _ e Q dW p 

p <P\n{d/Ry kB?o' {) 



where d is the lattice constant, R is the radius of the metallic wires, a is the 
static conductivity of the metal. The same results for e and oo p have been later 
obtained theoretically by Sarychev and Shalaev[8]. 

The San Diego group accepted the "plasma model" and considered the nega- 
tive e at uj < Up as one of the two crucial conditions for creation of the LHM. 
They have reported the first observation of the anomalous transmission and 
negative refraction in a compound system of split ring resonators (SRR's) and 
MPC[2,9]. According to Pendry et al. [10] such resonators create a negative 
bulk magnetic permeability due to the anomalous dispersion. 

In this paper we present an analytical theory of the EMW's propagating in 
the MPC without SRR's at very low frequencies. The theory is based upon the 
parameters /<C 1 and uR/c <C 1, where / is the volume fraction of a metal 
in the MPC. First we derive an exact dispersion equation for the s-polarized 
EMW in the system of infinite parallel thin straight wires ordered in a square 
lattice. The electric field of the wave is directed along the wires (z-axis), while 
the wave vector k is in the x-y plane. Assume that the total current in each 
wire is Io exp[—i(ut — k • r^)], where is a two-dimensional radius-vector of 
the wire in the x-y plane. The external solution for electric field E z of one 
wire, located at r; = 0, has a form 

K = (» , (3) 



where Hq = J — iN is the Hankel function which decays exponentially at 
Imw < 0. Neither J nor N has this important property. Here and below we 



2 



omit the time dependent factor. The solution is written in cylindrical coor- 
dinates z,p,(f> and it satisfies the boundary condition = (i/u))dE z / dp = 
I p /2np at p = R. 

The electric field created by all wires is 

E,(r) = ^e^Ee !k '< r -^?(^), (4) 



where pj = y (x — xj) 2 + (y — i/j) 2 , summation is over all sites of the square 
lattice and the sum is a periodic function of r. The dispersion equation fol- 
lows from the boundary condition [11] that relates the total electric field at 
the surface of any wire / to the total current through the wire E Z {R) = 
J exp (ik ■ ri)/a e fnR 2 , where a e f = 2a J\(kR) j 'kRJq(kR) , k = (l+i)/5, and 5 
is the skin depth. At small frequencies, when S > R one gets a e f ~ a. At high 
frequencies, when 5 < 0.1/2, one gets the Rayleigh formula a e f xs (1 +i)a8/R. 

Note, that the EMW exists mostly if the skin-effect in the wires is strong. 
Using Eqs.(l,2) one can show that Y/oo p = 5 2 /(R 2 \nd/R), where 5 is the skin 
depth at u = uo p (see also Ref. [12]). Thus, with the logarithmic accuracy one 
can say that r <^ u p if 5 <^ R. We show below that the exact solution has 
similar properties. 

Finally, the dispersion equation for u(k) = (x — i^)c/d has a form 

(X ~ il) E e^'+^W M = ^7, (5) 

l,m d(J eff 



where where zi m = (x — il)yl 2 + m2 + (R/d) 2 , I and m are integer numbers 
and the small term (R/d) 2 under the square root is important only when 
I — m — 0. Taking real and imaginary parts of Eq.(5) one gets two equations 
for x an d 7- 

Now we study the cutoff frequency for the EMW's uj , which is the solution 
of Eq.(5) at |k| = 0: 

(X-<7)E^ 2 W) = ^7- (6) 

The numerical results for the real part of the frequency are shown in Fig. 1. 
One can see that x is of the order of few units. The values of 7 are of the order 
of the right hand side of Eq.(6). Thus, x 3> 7 if Z^e/Ao^ ^> 1- 
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In addition to the numerical solution we propose an approximation valid at 
a very small /, when | ln/| ^> 1. We separate the term with / = m = in 
Eq.(6) and substitute the rest of the sum by the integral. Then 



(x - h) 




C 



4i 



4ce 



(7) 



where we assume that 7 C x- Here C is the Euler's constant. 



Figure 1 compares our result for the real part of the cutoff frequency given by 
Eq.(6) with the results given by Eqs.(2,7). In these calculations we assume that 
a e f is given by the Rayleigh formula so that the right hand side of Eqs. (6,7) has 
a form a(l —i)y/Xi where a = Ice^Rj j dab and S is taken at uo = c/d. One can 
see from Fig.l that approximation Eq.(7) is much better than approximation 
Eq.(2). Both approximations coincide at small / and they are accurate at 
extremely low values of / (~ 10~ 7 ) when the logarithmic term in Eqs. (2,7) 
is very large. The computational and experimental data of Refs.[5,9] are also 
shown at Fig. 1 and they are in a good agreement with Eq.(6). The result 
of Pendry group [6] for the cutoff frequency is 8.2 GHz, while Eq.(6) gives 
8.6 GHz. Experimental result of Bayindir et a/. [13] is 11.67 GHz while our 
Eq.(6) gives 11.8 GHz. Thus, we can make a conclusion that the San Diego 
group, group of Pendry and Soukoulis group discuss the same mode but at 
different values of parameters and that our analytical theory describes the 
same mode as well. 



At Fig. 2 we plot the imaginary part of the z-component of the electric field 
in the unit cell as calculated from Eq.(4) with k — 0. Real parts of the x- and 
y- components of the magnetic induction, are shown in Fig. 3. Real part of 
the E z and imaginary parts of B x , B y are much smaller than imaginary part 
of E z and real part of B x , B y correspondingly if 7 1. As one can see, the 
fields are strongly modulated inside the lattice cell and E z is close to zero near 
each wire so that absorption is small. This fact allows the low frequency mode 
propagate in the MPC under the condition fa e f/eou ^> 1. 

Now we find the component of electric permittivity e(u) = e zz , which describes 
the s-polarized extraordinary waves in the uniaxial crystal. It is defined by the 
relation e = e + iad/xc, where effective macroscopic conductivity a relates 
average current density J z to the average electric field E z by equation j z = oE z . 
To find a we introduce an external electric field £ z e~ lu)t . Using the boundary 
condition on a wire, which now has a form 

— 4d — h W+2 "^ (8) 
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one can find a relation between the current and the average field which gives 
both a and e. Finally one gets 



eo 



X 



-*x 



r«7 ^ TT {2) , .. N ce 
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fda ef f 



(9) 



The expression in the square brackets of Eq.(9) is the dispersion equation (6). 
One can see that Ree changes sign at u = u and becomes negative at u < u , 
where c^o is the root of the dispersion equation (6). 

To find a>(k) one should solve Eq.(5). For small |k| one can get an analytical 
result uj 2 = c 2 k 2 + uj 2 , which is isotropic in the x-y plane. 

Now we discuss the possibility of creation of the LHM using the negative e of 
the MPC. Suppose that the wires are embedded into a medium with the neg- 
ative magnetic permeability /i. One can see that in this case the propagation 
of any EMW is suppressed. Indeed, instead of Eq.(5) one gets equation 

J2 e *<*.*i+Mi>tf (^J x 2 + y 2 + R 2\ = 2 * 

j Vc v J J J \/j,\ujR 2 a e f 



where Cq = 1/y |/i|eo and K is the modified Bessel function. One can see that 
at |k| = all the terms on the left hand side of this equation are positive and 
real if Imw is small. Thus, if the right hand side is small, the equation cannot 
be satisfied. At small values of uo and |k| summation in the Eq.(lO) can be 
substituted by integration. Assuming that Reu ^> Imu one gets l + k 2 cl/u 2 = 
—ia e ff/uj€Q. This equation does not have real solutions for u(k). Thus, at 
negative /x there are no propagating modes at any frequency under the study. 

Now we compare this result with the theoretical idea[9,14] to obtain the LHM, 
where negative e is created by the system of wires and negative \x is created in 
some other way. This idea is based upon the assumption that the negative e 
at uj < ujq results from a "longitudinal plasma mode" . It is taken for granted 
that its frequency is independent of magnetic properties of the system, which 
is usually the case for plasmons. However, the mode discussed above is not a 
plasma mode (see also [15]). One can show that this mode has zero average 
value of the magnetic induction B over the unit cell. In this sense this is 
indeed a longitudinal mode. But the average value of the magnetic energy, 
which is proportional to B 2 , is not zero and it is large. The physics of this 
mode is substantially related to the magnetic energy. That is why negative \i 
completely destroys this mode. It destroys also the region of negative e. In fact 
this could be predicted from the observation that ujq ~ c/d becomes imaginary 
at negative /i. For example, one can see from Eqs.(l,2) that at /i < one gets 
uj 2 < and e > at all frequencies assuming that T is small. 
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Thus, we have shown that the simple explanation [14,9] of the negative re- 
fraction in the compound system of the MPC and SRR's, based upon the 
permittivity e of the MPC and the negative permeability \x of the SRR's does 
not work because negative \x blocks propagation of EMW's in the MPC. The 
propagation observed by the San Diego group might be a manifestation of the 
remarkable conclusion of Landau and Lifshitz (See Ref.fll] p. 268) that 
does not have physical meaning starting with some low frequency. Then, the 
explanation of the negative refraction in this particular system would be out- 
side the simple Veselago scenario (see Ref . [16] as an example) . In any case, to 
explain the negative refraction one should use a microscopic equation similar 
to Eq.(5) but with the SRR's included. 
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Fig. 1. The real part of the dimensionless cutoff frequency xo as a function of the 
volume fraction of metal /. Solid, dashed, and dotted lines represent solutions of 
Eq.(6), Eq.(7), Eq.(2) respectively. The experimental data of Ref.[9](-^") are shown 
together with numerical data of Ref.[5] (A, •, ■). On the main plot a = 0.024, 
d = 12.7 \xm for upper solid, upper dashed lines, and the point A; a = 0.078, 
d = \.11\xm for middle solid, middle dashed lines, and the point •; a = 0.246, 
d = 0.13/zm for lower solid, lower dashed lines, and the point ■. On the insert 
a = 3.4 • 10" 4 , d = 8.0mm. 




Fig. 2. Imaginary part of the electric field inside the lattice cell in the units of 
I^jiQc/d. The radius of the wires R = 0.8 mm, the lattice constant d = 8 mm, the 
static conductivity of the metal of the wires a = 5.88 • 10 7 (0 • m) _1 . 
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Fig. 3. Real part of the x-(a), and y-component(b) of the magnetic induction B 
inside the lattice cell in the units of Iona/d. The radius of the wires R = 0.8 mm, 
the lattice constant d = 8 mm, the static conductivity of the metal of the wires 
a = 5.88 • 10 7 (n ■ m) _1 . 
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